The longitudinal spin relaxation of 2d electrons in Si/SiGe quantum wells in a 

magnetic field 



O 

o 

(N 

o 
Q 

(N 



o 

-4-* 



o 
o 



> 



c3 



i 

C 

O 

o 



x 

S3 



Z. Wilamowski 1 ' 2 , W. Jantsch 1 
1 Institut fur Halbleiterphysik, Johannes Kepler Universitdt, A-4-040 Linz, Austria, and 
2 Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/4-6, PL 0668 Warsaw, Poland 

(Dated: February 1, 2008) 

The longitudinal spin relaxation time, Ti, in a Si/SiGe quantum well is determined from the 
saturation of the ESR signal. We find values of a few microseconds. Investigations of Ti as a function 
of Fermi energy, concentration of scattering centers and of the momentum scattering time, Tk, lead 
to the conclusion that for high electron mobility the spin relaxation is ruled by the Dyakonov- 
Perel (DP) mechanism while for low mobility the Elliott- Yaffet mechanism dominates. The DP 
relaxation is caused by Bychkov-Rashba coupling. Evaluation of the DP mechanism shows that 
Tj -1 for high electron mobility can be effectively reduced by an external magnetic field. The effect 
of the degenerate Fermi-Dirac statistics on the DP process is discussed. 



I. INTRODUCTION 

Looking for spin systems suitable for spintronics or 
quantum computing devices, the longitudinal spin relax- 
ation time, Ti, is of basic importance. T\ is ruled by 
spin-flip processes and it corresponds to the characteris- 
tic spin memory time. In this paper we investigate T\ 
in the high mobility 2d electron gas in a Si/SiGe quan- 
tum well, where electrons can be easily manipulated by 
illumination with light and by an electric field (l], |^] . We 
show that for this material system, which magnetically 
is one of the cleanest, T\ is of the order of microseconds 
whereas the time needed for a spin manipulation by a 
microwave magnetic field is by more than two orders of 
magnitude shorter. 

We also investigate the mechanism for spin relaxation. 
Analyzing the spin relaxation rate as a function of the 
momentum scattering rate, t^ 1 , allows to distinguish the 
Elliott- Yaffet (EY) mechanism || and the D'yakonov- 
Perel (DP) mechanism |fJ. The EY mechanism de- 
scribes the probability of a spin-flip in a momentum scat- 
tering event. This probability is ruled by spin-orbit cou- 
pling and the resulting admixture of a state with opposite 
spin projection fi 0: 



(ir 1 : 



EY = OlEYT, 



(1) 



The DP relaxation, in contrast, originates from a zero 
field spin splitting of the conduction band states ^, ||, |) . 
For a Si quantum well the zero field splitting is described 
by the Bychkov-Rashba (BR) term M §: 



Ti-BR = asi?(k x a) ■ e z 



(2) 



which was shown to exist also for single sided modulation 
doped Si quantum wells 

Here a stands for the vector of a Pauli spin-matrix of 
a conduction electron ||, k is the k- vector proportional 
to the electron momentum, e z is a unit vector perpen- 
dicular to the 2d layer and cxbr is the Rashba parameter 
that depends on the spin-orbit coupling and details of the 



interface |L(J. Momentum scattering causes also a time 
dependent modulation of the BR interaction. As a conse- 
quence, the probability for spin-flips becomes finite. For 
non-quantizing magnetic field, the DP mechanism j^, |(J is 
expected to be proportional to the momentum scattering 
time: 



(Tf 1 ) 



DP 



^BR T k 



(3) 



where the frequency Qbr is proportional to the k-vector 
and the BR parameter cxbr'- 



^br 2 = a BR k/2h 



(4) 



In this paper we present results obtained from con- 
duction electron spin resonance (CESR) spectroscopy. 
Simultaneous measurements of CESR, which allows to 
evaluate Ti, and cyclotron resonance (CR), which allows 
|l], |^] to estimate Tk, permit the evaluation of T\ as a 
function of T k . Such data are obtained from samples 
with different donor and electron concentrations. 



II. SAMPLES AND EXPERIMENTAL RESULTS 

Samples were grown by molecular beam epitaxy on 
1000 f2cm Si(001) substrates, which show complete car- 
rier freeze-out below 30 K. A 20 nm thick Si channel with 
tensile in- plane strain was deposited on a strain-relaxed 
Sio.75Geo.25 buffer layer, which consists of a 0.5 /Ltm thick 
Sio.75Geo.25 layer on top of a 2 /jm thick Sii-ajGe^ layer 
with compositional grading. The upper Sio.75Geo.25 bar- 
rier was modulation doped with a 12.5 nm thick, nomi- 
nally undoped spacer layer, and capped with 5 nm of Si. 
Three modulation doped Si/SiGe structures with differ- 
ent donor concentrations were examined. The electron 
concentration was changed by the light illumination. 

All measurements were performed with a the standard 
X-band ESR spectrometer, at a microwave frequency 9.4 
GHz. The sample were situated in the center of the rect- 
angular TM201 cavity, at the maximum of the magnetic 
component of microwave field (which is perpendicular to 
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FIG. 1: (a) Longitudinal spin relaxation time of the 2d 
electron gas in a Si/SiGe quantum well as a function of the 
sheet electron concentration, n s . (b) Concentration depen- 
dence of the momentum scattering rate, t^T 1 , as evaluated 
from the cyclotron resonance linewidth. Different symbols 
stand for samples with different donor concentrations in the 
doping layer. 



the applied magnetic field). The sample layer was ori- 
ented to be perpendicular to the applied magnetic field 
(and to electric component of microwave field). 

The spin resonance has an exceedingly narrow 
linewidth in the range 3 -j- 10//T. In spite of the fact that 
the sample was situated in the minimum of the electric 
microwave field and perpendicular to it, the strong ab- 
sorption due to CR was well observed allowing for to 
monitor the carrier density from the integral absorp- 
tion and the momentum scattering rate from the CR 
linewidth j|. 

In Fig. 1 the longitudinal spin relaxation time, T%, 
and the momentum relaxation rate, tZ\ are plotted as 
a function of the electron concentration, n s . The param- 
eter T\ has been evaluated from the saturation of the 
ESR signal amplitude and the ESR line broadening at 
high microwave power Estimating the quality 

factor of the loaded cavity from the resonance dip width, 
we obtain the amplitude of the magnetic component of 
the microwave field of 1.1 G at a microwave power of 200 
mW. The data for different samples are marked by dif- 
ferent symbols. The results for the spin relaxation time, 
Ti, vary in the range of 1 to 5 /is. For different samples 
T\ is different and it depends on the electron concentra- 
tion. The momentum scattering rate varies with n s by 
an order of magnitude. The increase of the momentum 
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FIG. 2: (a) Spin relaxation rate, Tj -1 , as a function of mo- 
mentum scattering rate, . It shows that for the strong 
scattering the EY process dominates the spin relaxation. The 
spin relaxation normalized by the electron concentration, n a , 
is shown in Fig (b). It shows that the spin relaxation for the 
high electron mobility, t^ 1 < 3 • 10 11 s _1 , is dominated by 
the DP mechanism. Solid line: DP relaxation rate for a BR 
parameter of obr = 1.1 • 10 11 eVcm || and degenerate statis- 
tics (Eqs.(7-8)). Dash-dotted line: non-degenerate statistics 
(Eqs. (6-7)). Dashed straight line: from Eq. (2) which de- 
scribes the DP relaxation for non-degenerate statistics in the 
absence of an external magnetic field. 



scattering is related to the screening breakdown and an 
increase of the potential fluctuations at low Fermi energy 
Q . Samples with a higher doping level show also a higher 

T k ■ 

The dependence of T± on n s is governed by the complex 
dependence of the relaxation rate on the Fermi k-vector 
and of the dependence of t^T 1 on the electron concen- 
tration. In order to follow the dependence of the spin 



relaxation rate 



T -i 



on the momentum relaxation, ty. 



our data are plotted in Figs. 2a and 2b in two different 
ways. In Fig. 2a the spin relaxation rate, Tj -1 , is given as 
a function of the momentum scattering rate, t^T 1 . In Fig. 
2b, the spin relaxation rate, T^ 1 , is normalized by the 
electron concentration, n s . This normalization allows to 
account for the dependence of T^ 1 on the BR parameter, 
a br, and to study the dependence of the DP rate T^ 1 
on TjT 1 . 
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III. THE SPIN RELAXATION CAUSED BY 
THE ELLIOTT- YAFET MECHANISM 

Comparison of the data in the two figures demonstrates 
the existence of two different ranges with different spin 
relaxation behavior. For high scattering rate, r^ 1 > 3 • 
I 11 s^ 1 , the spin relaxation is simply proportional to the 
momentum scattering rate indicating the EY process as 
the dominating one. The EY coefficient is independent 
both of the electron concentration and the doping level 
within the experimental error. The solid line in Fig. 2a 
corresponds to a^y = 1.0 • f0~ 6 . 

For low momentum scattering rate, TjT 1 < 3 ■ I0 11 s — 1 , 
the spin relaxation rate is bigger than expected from the 
EY mechanism. Moreover, T{~ depends on the electron 
concentration in the high electron mobility range. On the 
other hand, the normalized spin relaxation (see Fig. 2b) 
is characterized by a systematic dependence, common for 
all investigated samples. As we argue below the observed 
dependence in the high mobility range is well described 
by the DP mechanism. 



IV. THE D'YAKONOV-PEREL MECHANISM 
OF THE RELAXATION 

The prediction of the DP scattering rate, as described 
by Eqs. (3) and (4) is marked by the dashed line in 
Fig. 2. For the BR parameter we took the value 
OiBR = 0.55- I0~ 12 eV cm evaluated earlier from the anal- 
ysis of the linewidth and g-factor anisotropy in the same 
samples J^]. No correlation between the dashed line and 
the experimental data is recognizable. Eq. (3) stands, 
however, for the case of a weak external magnetic field, 
when the momentum scattering rate is much smaller as 
compared to the Zeeman frequency, and it does not con- 
sider cyclotron motion and Landau quantization. In an 
external magnetic field, because of the cyclotron motion, 
the electron velocity changes its direction all the time. 
The time correlation function of the k-vector, and conse- 
quently the correlation function of the effective BR field 
seen by an electron, is describes by: 



(k • k(t)) ex (n BR • n BB (t)> = n 2 BR e l - 



(5) 



The corresponding probabilities of spin-up and down 
flips are obtained by the Fourier component of the cor- 
relation function (Eq.(5)) at the Zeeman frequency lo : 



W± = 



n 2 



Tk 



2 1 



((Jr. ± UJ, 



(0) 



The longitudinal spin relaxation time for a single elec- 
tron is equal to: 



(T^) DP = W + + W^ 



(7) 



Eqs. (6) and (7) describe the DP relaxation in an ex- 
ternal magnetic field. For a short momentum relaxation 



time, Eq.(7) becomes equivalent to Eq.(3). But for quan- 
tizing magnetic field, where lo c t^ > f , the DP relaxation 
rate is expected to be reduced by the denominator in 
Eq. (6). The dependence corresponding to Eqs. (6-7) is 
shown in Fig. 2 by the dash-dotted line. The reduction 
of the spin relaxation caused by the external magnetic 
field is well visible. Moreover, for low scattering rate, 
<-L>cT~k 3> f , the DP relaxation rate is expected to be pro- 
portional to the momentum scattering. 

For an electron gas, in which the final states to which 
electron can be scattered, are partially occupied the eval- 
uation of the mean spin relaxation rate of the whole 
electron system requires thermodynamic averaging. The 
scattering probability and the momentum relaxation rate 
depend on energy. These quantities are proportional to 
the population of empty states: r / 7 1 (e) = t^ [1 — fpo (e)] 
where t^ d is the momentum relaxation rate (as used in 
the Boltzmann equation approach) and Jfd{£) is the 
Fermi-Dirac distribution function. For moderate mag- 
netic field, where Huj c < k B T < Ep the dependence of 
the BR frequency on energy can be neglected and the 
mean value of the transition probability, (W± } , weighted 
by the derivative of the Fermi-Dirac distribution func- 
tion, is described by: 



n 2 



Tfe(e) 



f + (lu c ±uj ) 2 t 2 (s) 



de (8) 



The solid line in Fig. 2b corresponds to the DP scat- 
tering rate as described by Eqs. (7) and (8). For the BR 
parameter again a value of 0.55 • I0~ 12 eV cm has been 
taken. The theoretical curve, without any other fitting 
parameter, fits well to the experimental points for high 
electron mobility. For the highest mobility the effect of 
a moderate magnetic field (B=0.34 T) is the reduction 
of the DP mechanism by about two orders of magnitude. 
In the limit of very high electron mobility, the DP re- 
laxation rate for degenerate statistics (Eq. (7-8)) is by a 
factor 2 smaller as compared to the non-degenerate case 
described by Eqs. (6-7) 

For the high momentum scattering rate the DP relax- 
ation rate is expected to tend to the solution for weak 
magnetic field. But for degenerate statistics (solid line), 
where the final states are partially occupied, the spin re- 
laxation rate is by a factor 2 bigger as compared to the 
non-degenerate case (dashed and dash-dotted lines) for a 
given momentum relaxation rate. 



V. CONCLUSIONS 

In conclusion, we have shown that: 

• the DP mechanism dominates for high mobility 
structures but the quantization due to the applied 
magnetic field leads to a considerable reduction of 
the DP relaxation rate. In Eq. (3) a reduction fac- 
tor of about I + t 2 uj 2 must be introduced, where 
uj c is the cyclotron frequency (compare Eq. (3) and 
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Eqs. (6-7)). As a consequence, for weak momen- 
tum scattering the reduced DP spin relaxation rate 
is proportional to r^T 1 , in contrast to Eq. (3). 

• the value of the BR parameter, asn, as determined 
from T\ turns out to be the same within the exper- 
imental accuracy as previously evaluated from the 
anisotropy of the CESR linewidth (dephasing time, 
T 2 ) and the g-factor [|. 

• for low mobility samples the spin relaxation is dom- 
inated by the EY mechanism. We find an EY co- 
efficient of: oley = 1-0 • 10~ 6 , which is common 
for all samples and, for the investigated range of 



parameters (the Fermi energy does not exceed 2.5 
meV), c*ey does not depend on the electron con- 
centration. 
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